Current partition at topological zero-line intersections 
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An intersection between one-dimensional chiral acts as a topological current splitter. We find 
that the splitting of a chiral zero-line mode obeys very simple, yet highly counterintuitive, partition 
laws which relate current paths to the geometry of the intersection. Our results have far reaching 
implications for device proposals based on chiral zero-line transport in the design of electron beam 
splitters and interferometers, and for understanding transport properties in systems where multiple 
topological domains lead to a statistical network of chiral channels. 



A massive chiral two-dimensional electron gas 
(C2DEG) has a valley Hall conductivity that has the 
same sign as its mass. The valley Hall effect leads to 
conducting edge states and also, when the mass param- 
eter varies spatially, to conducting states localized along 
mass zero-lines. — ^— Provided that inter- valley scattering 
is weak, zero-line state properties are closely analogousi"^ 
to edge state properties of quantum spin-Hall insula- 
tors and include both chiral propagation and suppressed 
backscattering."'^ Metallic zero-line modes (ZLMs), or 
topological ID kink states, provide a two dimensional 
realization of Dirac zero energy modes^i^ and their ex- 
istence has been proposed in a wide variety of systems 
including graphene mono and bilayerSf^"— i^"— topologi- 
cal insulators with lattice dislocations?^ boron nitride 
crystals with grain boundaries/- superfluid '^He,»ii and 
photonic crystals ji^ii^ In the present Letter we examine 
current partition properties at zero-lines intersections)^ 
which are expected to be ubiquitous in systems in which 
the mass term results from a disorder potential or from 
spontaneous symmetry breaking. 

ZLMs in C2DEGs are centered on zero-lines of the 
massi"— i^ii, i.e. on lines along which the mass changes 
sign as illustrated in Figure [T^. A mass term leading to 
a valley Hall effect ^'''^^ can be produced by a sublattice 
staggered external potential in single layer graphene,-^ 
and more practically by a gate controlled interlayer po- 
tential difference in Bernal bilayer and ABC stacked mul- 
tilayer graphene. ^r^i^ Mass terms can also be generated 
by spin-orbit couplingi^"— and by electron-electron inter- 
actionsji^Ti^ In this last case ZLMs^ appear naturally 
at domain walls separating regions with different local 
anomalous, spin, or valley Hall conductivities. 

Chiral propagation implies that ZLMs can travel only 
in the direction which places negative masses either on 
their left, or depending on valley, on their right. It fol- 
lows, as illustrated in Figure [T]:, that there is no for- 
ward propagation at a zero-lines intersection; a propa- 
gating mode is split between a portion that turns clock- 
wise and a portion that turns counterclockwise. These 
unusual transport properties are potentially valuable for 
new types of electronic devices. We have therefore carried 
out quantum transport calculations for an explicit model 
of intersecting ZLMs in order to discover rules for current 



partitioning at such a ZLM splitter. The system we study 
is a TT-band tight-binding model for single-layer graphene 
with a position-dependent sublattice-staggered potential 
constructed to form intersecting zero lines which enable 
propagation to four ZLMs labeled left (L), right (R), 
up (U) and down (D) in Figure [TJj. For simplicity, we 
consider the case where the U and D ZLMs are fixed along 
the vertical direction, and we define the angles between 
R and D ZLMs to be a and between L and U ZLMs to 
be p. The blue and orange lines in Figure [TJd indicate 
the allowed chiral propagation paths. 

The numerical results reported on below are for a 
TT-orbital tight-binding Hamiltonian with nearest neigh- 
bor hopping and a sublattice-staggered potential: H = 

C^CJ+t7AEi,AC^C^ + C^BE^,BC^C^■ Here cj {Ci) 

is a creation (annihilation) operator for an electron at site 
i, and t = 2.6 eV is the nearest neighbor hopping ampli- 
tude. For a sublattice staggered potential the A and B 
sublattice energies are opposite, i.e. Ua = —Ub — AC/q, 
where Uq measures the potential strength and A="±" 
determines the sign of the valley Hall effect in each quad- 
rant. In all our simulations, the potential amplitude was 
chosen to be Uo/t = 0.05. The ZLMs appear confined 
along zero-lines where the mass becomes zero and the 
resulting wave-function tails spread into the bulk with a 
depth proportional to the inverse of the massi^i^. 

Although our study has been carried out in single 
layer graphene for computational convenience, we ex- 
pect that similar conclusions apply to ZLMs in other 
systems, and in particular in bilayer graphene where the 
spatially patterned mass term can in principle be gener- 
ated externally using gates. Our transport calculations 
are based on the Landauer-Biittiker formalism^ and re- 
cursively constructed Green functions. ^^'^^ The conduc- 
tance from lead q to lead p is numerically evaluated 
from Gpq = (eV/i)Tr[rpG'TqG°], where e is the elec- 
tron charge, h is the Planck's constant, C'" are the re- 
tarded and advanced Green's functions, and Tp is a 
linewidth function describing the coupling between lead 
p and the central region. The propagation of ZLMs, 
or kink states, incoming from lead p is effectively illus- 
trated by plotting a map of its contribution to the lo- 
cal density of states (LDOS) at an energy e in the gap: 
Pj,(r, e) = l/27r[G''rpG°]rr- Here r is the real space co- 
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FIG. 1: Current partition at zero-line intersections: a. Staggered sub lattice-potential, and hence valley Hall conductivity, 
spatial distribution pattern which defines the zero-line paths. Blue dots indicate positive site energies and green dots negative 
site energies, b. Schematic of four-terminal graphene samples with different staggered potential distributions. The Left (L), 
right (R), up (U), and down (D) leads are extend indefinitely from the plotted central scattering region. The U and D ZLMs 
are fixed along the vertical direction and the angles a and /3 specify the R and L ZLM directions. The solid and dashed lines 
in black denote zero lines. The blue and orange lines represent allowed chiral propagation paths, c. ZLM LDOS distribution 
for modes incident from lead R for a = 90°, 60°, 120° at fixed /? = 90°. 



ordinate. 

The central scattering region in our calculations is rect- 
angular with size Ux = 94 and Uy = 432 as explained fully 
in the Supplementary Information. The valley label of a 
state is of course not a good quantum number for ZLMs 
and valleys are most strongly mixed when their wave vec- 
tor projections in the propagation direction are identical. 
For energies inside the gap this coincidence happens only 
for propagation in the armchair directionii^ii^iiS. How- 
ever, numerical calculations have shown a remarkable 
absence of back-scattering at sharp turns in the zero- 
line or at the encounter of a ZLM splitter— except in a 
narrow energy range very close to the avoided crossing 
gap centered on e/t = 0.00 between modes with opposite 
propagation directions. For the results shown below we 
have chosen e/t = 0.01 to avoid this energy range; the 
chirality of the ZLM modes is then very well defined. 

In a four terminal ZLM splitter device [see Figures [TJi], 
there are in total twelve distinct inter-terminal conduc- 
tance values. The number of independent conductances 
is reduced to six in time-reversal symmetric systems 
since Gpq = Gqp. For chiral transport forward scat- 
tering and back scattering are absent at a ZLM inter- 
section, reducing the number of independent parame- 
ters further. The current conservation then implies that 
Gpr + Gqr = Gq = c^/ft. for any value of r and p, q the 
labels of the two neighboring leads. It follows that 

Glu = Grd & Gru — Gld, (1) 

and that Gru + Grd = e^/h, leaving only one indepen- 
dent parameter for the entire four terminal systems. In 
a ZLM splitter with zero backscattering and perfect chi- 
ral current filtering transport is completely characterized 
by specifying how incoming current at an intersection is 



partitioned between clockwise and counterclockwise ro- 
tation outgoing directions. The partition law must be the 
same for all incoming channels. In the following, we fo- 
cus on the conductances Gur and Gdr corresponding to 
the currents incoming from lead R. The above relations 
were numerically verified for a ZLM current splitter with 
a = (3 = 90° in Ref. [3| and we have now numerically 
verified that they are true for arbitrary values of the lead 
angles a and (3. (See the Supplementary Information for 
further details.) 

For a ZLM splitter with a fixed vertical pair of leads U 
and D and rotatable R and L leads as shown in Figure[T}D, 
we have numerically discovered a rather simple law which 
describes the dependence of the current partition on the 
angles a and (3 with surprising accuracy: (See Figure [2] 
for a summary of the data which supports this law and 
the Supplementary Information for further details.) 

Gur = ^[l-sin(a + /3)], 

Gdr = ^ [I + sm{a + 13)], (2) 

for 90° <a + l3 < 270°. Outside of this angle range the 
current follows the path with the larger rotation angle: 
Gdr = GuL = Gq and Gur — Gdl — when 0° < a + 
^ < 90°, or Gur = Gdl = Go and Gdr = Gul = when 
270° <a + /3 < 360°. For the special case oi (3 = 90° the 
current partition law simplifies to Gur = Gosin^(Q;/2) 
and Gdr — Go cos^(q;/2). And for the special case of 
a — /3, corresponding to zero lines that are straight at the 
intersection point and therefore to a mass pattern defined 
by a smooth external potential, Gur — Gq sin^(a — 45°) 
and Gdr = Go cos'^{a - 45°). 

Gontrary to intuition, as illustrated in Figure [Ij; for 
a current incoming from lead R, the current partition 
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FIG. 2: Current partition law: Current as a function of 
angles a and /3 for the zero-line intersection defined in Fig. 
Ic. The current partition follows a simple relation given by 
Eq. [2] whenever 90° < Q + /3 < 270° and completely follows 
the sharper rotation path outside this range. The solid and 
dashed lines are the fitting functions defined in Equation 



law implies that more current always follows the path re- 
quiring a larger rotation angle. This makes sense since 
this pair of incoming and outgoing modes interact over a 
longer distance, providing more opportunity for interac- 
tions which lead to inter-mode scattering (The zero-line 
mode decay length is ^ mvD/fi where m is the local 
mass and is the Dirac point velocity for m = 0.—). It 
is noteworthy that the current partition depends only on 
the combination a + /3 and indeed it is surprising that 
the angle (3, which specifies the orientation of a lead that 
does not carry any current, influences the partition law. 
This law evidently results, however, from the interfer- 
ence between ZLMs close to the ZLM intersection point, 
allowing inactive leads to play a role in the scattering. 
Note that the completely unbalanced current partition 
favored by a very small value of a (a sharp turn from R 
to D) is mitigated if (3 (the L to U turn angle) is larger 
than 90°. 

We emphasize we have not derived the form which fits 
our current partition numerical results analytically and 
it is likely only approximate. It applies accurately only 
at energies close to the middle of the bulk energy gaps. 
The expression can be rationalized by the following ar- 
gument. Let us consider a ZLM splitter, where the chiral 
current filtering rule restricts current incoming from R 
or L to scatter into the U and D leads. By denoting the 
asymptotic amplitudes of the ZLM at the leads as An, 
Al, Ajj, and Ajj, current conservation implies that 
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Equation ([3]) relates the amplitudes of incoming and 
outgoing waves via a general SU(2) unitary transforma- 
tion matrix with parameters u, v, and r. When we set 
Al to be zero, then for any Ar, = |Aij|^sin^T, 
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FIG. 3: Influence of disorder on the current partition 
law: Averaged conductances (Gdr) (in panel a) and (Gur) 
(in panel b) from lead R to leads D and U as a function of 
disorder strength W at fixed /3 = 90° . Four representative an- 
gles of incidence a — 30°, 60°, 90°, and 120° are considered. 
Over 50 samples are collected at each point. The error bar is 
used to indicate the strength of fluctuations as a function of 
disorder realization. 



and l^cP — \Aji\'^ cos^ T. The scattering phases u and 
V are irrelevant for the scattering probability and will 
henceforth be dropped. The conductances from R can 
thus be written as Gur = Go sin^ r(a, (3) and Gdr — 



Go cos T {a, 13), where Go 



^/h. Equipartition for 



{a, 13) = {a, 180° - a) implies that T{a, 180° - a) = 45°. 
For the current coming from R into D lead, we get 
from symmetry considerations the relation GDR(a,/3) = 
Gdr(/5, ct). If we assume that t is a smooth function of a 
and /3 and additionally require that the conductances are 
invariant under {a, (3) ^ {a + N x 360°, P + M x 360°) 
we can conclude that T{a, (3) — c{a + f3)/2 — 45°, where 
c = ±1. Then we impose a second condition for the cur- 
rent partition saturation when a+/3 = 90°, for instance in 
the limit when /3 = 90° and a = 0° where R and D ZLMs 
merge together, that further restricts c = -1-1. This ex- 
pression reproduces Equation ^ for Gdr, and Gur, out- 
side the saturation region which occurs ioi a + /3 < 90° 
and a + 13 > 270°. 

Now we examine the robustness of our results in the 
presence of disorder. Long-range disorder is not effec- 
tive in producing inter valley scatteringi so we focus 
on the potentially more important short-range disor- 
der which we model as a random potential at each lat- 
tice site. Specifically we add a term Hdis — Si'^icjci 
with uji being uniformly distributed in the interval 
of [—W/2, +W/2], where W characterizes the disorder 
strength. In our simulations we considered over 50 real- 
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izations of the random disorder potential for each value 
of the strength. Figure [3] illustrates our results for the 
average conductances (Gur) and (Gdr) as a function 
of disorder strength W. We see that the current parti- 
tion law remains accurate up to disorder strengths larger 
than the bulk band gap A/t = 0.1. When the disorder 
strength is further increased, the averaged conductances 
(Gur) and (Gdr) are slightly reduced and disorder fluc- 
tuations grow. For example, when the disorder strength 
reaches W/t = 1.0, 10 times larger than the bulk band 
gap, the averaged conductances still retain over 80% of 
their original values. The lost current takes advantage of 
the intervalley scattering to access the outgoing modes 
of the L lead or to backscatter in the R lead. All these 
findings strongly indicate that the current partition law 
is very robust to disorder, suggesting that ideal zero-line 
transport properties can be approximated in real devices. 

Since there are presently no practical techniques for 
imposing staggered sublattice potentials in single layer 
graphene, other closely related systems may ultimately 
be of greater experimental interest. In Bernal stacked 
bilayer graphene, for example, a ZLM splitter can be re- 
alized by using gates to achieve perpendicular electric 
fields which vary in sign spatially. Another possibility 
is 2D honeycomb photonic crystals, in which the Dirac 
points have been experimentally observed and sublattice 
staggered potentials can be realized by choosing different 
diameters for the cylinders which form the structure or 
by varying the dielectric material used. 

In summary, when intervalley scattering can be ne- 
glected, transport along the zero lines of a sublattice- 
staggered potential in graphene is chiral, requiring travel 
in a direction which keeps positive masses on either the 



left or the right, depending on valley. We have used the 
Landauer-Buttiker formula and recursively constructed 
Green's functions to examine how chiral currents are par- 
titioned between available outgoing leads at a ZLM inter- 
section. We find that at energies near the middle of the 
bulk gap our numerical results for the dependence of cur- 
rent on ZLM geometry are accurately described by a sim- 
ple partition law specified in Eq. ([2]), and that the infiu- 
ence of disorder on this law is weak. The helicity of ZLM 
provides a new mechanism for allowing or blocking cur- 
rents and may find applications in alternative designs for 
nanoelectronic devices or in enabling electron quantum 
interferometry^^ in a new setting. We have explored, for 
the first time to our knowledge, the geometry-dependent 
current partition laws at the intersection of two zero- 
lines. It will be interesting to extend our present studies 
to more general parameter spaces and to look for similar- 
ities and differences with respect to other systems with 
chiral ID transport channels including photonic crystals, 
quantum anomalous Hall and quantum Hall effect sys- 
tems, and chiral superconductors. 
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